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An approach to nonlinear dynamical modeling of interaction between 
conscious and automatic processes in the brain is proposed. Illustration of 
this approach on the nonlinear equation for the current density in the cortex is 
presented. Nonlinearity is determined by the sigmoidal firing rate of neurons. 
The current density is considered as complex field with real and imaginary 
parts representing automatic and conscious processes. The interaction is due 
to the nonlinearity of the system. From this approach it follows that the 
nonlinear dynamics for pure imaginary process (a dream) is quite different 
from the dynamics of automatic process. 



Ubiquitous phenomena are often the most difficult to explain. During the 
last two centuries the problem of consciousness (C) was considered of limits 
for any scientific explanation. However, recently a number of prominent 
scientists took this problem seriously [1-3]. There are regular conferences 
on C, which attract hundreds of neuroscientists, physicists, mathematicians 
and philosophers. The possible advantages of developing a theory of C are 
enormous: from medicine to sophisticated robots to pacification of religious 
confiicts. 

From the physical-mathematical point of view it is desired to catch C into 
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some sort of equations. It is common knowledge that C is somehow connected 
with the electrochemical activity in the brain. So, it seems logical to start 
with equations for these processes. The brain activity revealed the regime of 
self-similarity [4] , which is typical for systems with strong interaction of many 
degrees of freedom. Corresponding equations can be formally written and are 
very complicated. However, at this stage the precise form of the equations is 
not critical. One can use various simplified models. The important question 
is how to connect these equations with C ? 

The C-processes are subjective and, as far as we know, they can not be 
measured directly by the objective methods, which are used for measuring 
electrochemical (automatic) processes. At the same time, there are reasons 
to believe that C-processes can interact with the automatic {A) processes. 
We need equations for A-fields and C-fields, which interact despite the fact 
that C-fields have a different nature and can not be measured by the same 
methods as ^-fields. In this Letter we suggest that processes in the brain can 
be described by generalized (C) complex fields, which have real A-component 
and imaginary C-component. These components interact due to the nonlin- 
earity of the equations. Such approach seems to be the simplest. If necessary, 
more general constructions can be considered in the future and they include 
hypercomplex components (quaternions), additional space (branes [5]), ad- 
ditional time [6] and distributed sources [7] . 

As first approximation, we assume that equations for the C-fields have 
the same form as equations for the A-fields. Let us illustrate this with a 
simple example. Consider model equation for the average (spatially uniform) 
current density a{t) perpendicular to the cortical surface: 

^ + /ca = /(a + cr) (1) 

Here k is the relaxation coefficient, a{t) is the average sensory input and / 
represents the sigmoidal firing rate of neurons (for example, f{a) — tanh(Q;)). 
For the case of spatially nonuniform Q;(t, x) and (7(i, x) we can use a more 
general equation, which include typical propagation velocity of signals v. 
Time differentiation of (1), simple manipulation and addition a term with 
the two-dimensional spatial Laplacian A gives: 

ol da d 

-g^ + {k + m)— + {km - v'A)^ = (m + ^)/(tt + (^) (2) 

where m is an arbitrary parameter. This type of equations are used for 
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interpretation of EEG and MEG spatial patterns (see recent paper [8] and 
references therein). In this context we have parameters : /c ~ m ~ v/i, where 
/ is the connectivity scale. 

Returning to (1), we now introduce the G-field g = a + iip, where ip 
represents the C-effect. Substitution of g instead of a into (1) gives system 
of two equations: 

da 

— + ka^Re{f{a + i'i/j + a)} (3) 

L/t' 

dib 

+ki^ ^ Im{f{a + ii^ + a)} (4) 

These equations are coupled because f{a) is nonhnear. Thus, we got sim- 
ple scheme for the A — C interaction. The same scheme can be apphed 
to equation (2) and to any nonlinear model equation. Note, that so-called 
extra-sensory effects (if they exist) can be included in this approach by as- 
suming that a has an imaginary part. When a and a are relatively small (in 
a dream), asymptotically (4) gives closed equation for ^0. The nonhnear term 
I'm{f{iip)} is quite different from the nonlinear term f{a) for the 74-process. 
For example /(a) = tanh(Q;) gives Im{f{iil))} = tsin{i/j). This may explain 
why dreams have bizarre dynamics, not only in content, but also in intensity. 

A lot of questions can be asked about the proposed approach. For in- 
stance, how C-fields are related to our emotions, thoughts, images etc.? Note, 
that we are modeling only an "active" part of C-field, which interact with 
certain A-ficld. It is important to design experiments in which a part of 
C-field is affected without substantial direct effect on corresponding A-field. 
C-modeling, which is related to electromagnetic activity in the brain, can be 
tested by using multi-channel EEG and MEG. More general C-modeling can 
include connection between electromagnetic activity and blood flow in the 
brain. In this case, brain imaging (particularly, MRI) can be used. There 
are some technical issues involved in such testing. But it seems doable. 

Another important aspect of C-modeling is its connection with the first 
principles. Introduction of imaginary fields (along with real fields) produces 
generahzations, which potentially can be used for treating some "ghosts" [5,9] 
in theoretical physics. Experimentally, it is intriguing to search for materials 
(besides the brain tissue), which can support an imaginary field. 
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